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Abstract 

To  quantify  the  effects  of  interactions  between  various  microstructure  attributes  on  fatigue  life  in  the 
high  cycle  fatigue  (HCF)  regime,  we  have  proposed  a  new  microstructure-sensitive  extreme  value 
statistical  framework.  This  framework  couples  the  extreme  value  distributions  of  certain  fatigue  indicator 
parameters  (FIPs)  or  response  functions  to  the  correlated  microstructure  attributes  that  exist  at  the  extreme 
value  locations  of  these  FIPs.  We  demonstrate  the  application  of  this  statistical  framework  to  investigate 
the  microstructure-sensitive  fatigue  response  of  the  PM  Ni-base  superalloy  IN  100.  To  accomplish  this  task, 
we  construct  statistical  volume  elements  (SVEs)  used  to  compute  the  local  response  for  200  instantiations 
of  IN  100.  These  SVEs  are  constructed  and  simulated  via  the  finite  element  method  with  crystal  plasticity 
constitutive  relations.  The  results  of  the  simulations  are  used  to  explore  extreme  value  statistics  of  the  FIPs 
for  these  microstructures.  The  extreme  value  distributions  of  the  Fatemi-Socie  FIP  are  fit  with  high 
confidence  by  the  Gumbel  distribution  and  are  defined  in  a  representative  nature  with  as  few  as  25 
simulated  microstructure  instantiations  (i.e.,  SVEs).  The  extreme  value  marked  correlation  functions  of  the 
apparent  Schmid  factor  based  on  the  geometry  of  the  slip  systems  relative  to  the  loading  direction  indicate 
that  cube  slip  may  be  important  to  fatigue  crack  formation  in  this  material  system.  This  supports  previous 
experimental  observations  of  fatigue  crack  formation  and  microstructurally  small  fatigue  crack  growth 
along  cube  planes  in  IN  100  in  grains  that  are  unfavorably  oriented  for  octahedral  slip  at  elevated 
temperatures. 

Keywords:  High  Cycle  Fatigue,  Ni-base  superalloys,  extreme  value  statistics,  probabilistic  fatigue,  and 
microstructure-sensitive  design. 
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1.  Introduction 

Scatter  in  the  High  Cycle  Fatigue  (HCF)  life  of  specimens  or  components  depends  on  the  extreme 
value  probabilities  of  having  existing  hot  spots  or  regions  with  increased  local  driving  forces  for  fatigue 
damage  formation  (i.e.,  fatigue  crack  formation  and  microstmcturally  small  crack  propagation). 
Specifically,  the  probability  of  fatigue  damage  formation  in  a  particular  volume  of  material  is  established 
by  the  extreme  value  (i.e.,  rare  event)  probability  of  a  particular  existing  combination  of  microstructure 
attributes  that  couple  with  the  applied  stress  state  such  that  fatigue  cracks  form  and  propagate  fatigue 
cracks.  Coupling  of  microstructure  attributes  with  loading  conditions  and  the  resulting  fatigue  response  is 
the  main  source  of  both  scatter  and  size  effects  in  fatigue  in  the  absence  of  other  random  environmental 
factors  (e.g.,  temperature,  atmosphere).  As  such,  we  commonly  observe  scatter  in  the  fatigue  response 
between  multiple  material  volumes  (or  components)  even  when  they  are  fabricated  from  the  same  batch  of 
processed  material  and  tested  in  nearly  identical  environments.  Although  environmental  effects  can  also 
contribute  to  the  probability  of  surface  to  subsurface  transition  of  crack  formation,  they  are  not  considered 
in  this  particular  work.  The  main  emphasis  of  this  work  is  on  developing  a  consistent  methodology  to 
characterize  correlations  of  microstructure  attributes  that  exist  with  high  probability  relative  to  the  nominal 
microstructure  in  regions  where  the  fatigue  driving  forces  are  maximum  (i.e.  extreme  value). 

The  dependence  of  fatigue  damage  formation  on  various  microstructure  attributes  has  been 
investigated  extensively;  several  reviews  have  been  published  on  the  subject  [1,  2].  For  example,  second 
phase  inclusions  or  pores,  large  grains  that  are  favorably  oriented  for  slip,  or  grains  adjacent  to  grains  that 
are  unfavorably  oriented  for  slip  have  each  been  linked  to  fatigue  damage  formation  in  several  metallic 
polycrystalline  material  systems.  Identifying  the  underlying  microstructure  attributes  that  drive  fatigue 
damage  formation  in  material  systems  with  multiple  phases,  complex  processing  histories,  etc.,  however,  is 
often  complicated  by  the  fact  that  multiple  interacting  microstructure  features  couple  with  the  imposed 
cyclic  plastic  deformation  and  stress  state  to  increase  the  local  driving  forces  for  fatigue  damage  formation. 
In  most  cases,  this  dependence  cannot  merely  be  deduced  from  direct  quantitative  image  analysis  of  various 
microstructure  attributes.  Assessment  of  the  coupling  of  microstructure  attributes  with  the  driving  forces 
for  fatigue  damage  formation  requires  a  combination  of  experiments  and  computational  simulation,  along 
with  a  connective  framework  based  on  extreme  value  statistics. 
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In  general,  fatigue  damage  formation  in  polycrystalline  metallic  materials  is  primarily  driven  by 
irreversible  slip  on  the  scale  of  the  microstracture.  In  addition,  interacting  microstructure  attributes  (e.g., 
grains,  phases,  inclusions,  and  voids)  can  increase  local  slip  and  associated  driving  forces  for  fatigue 
damage  at  the  microstructure  scale.  In  HCF,  where  stress  amplitudes  remain  below  the  macroscopic  flow 
stress  and  cyclic  plasticity  is  quite  heterogeneous,  fatigue  lives  are  dominated  by  fatigue  damage  formation 
rather  than  by  physically  small  or  large  crack  propagation  (cf.  [3-6]).  Accordingly,  the  emphasis  is  placed 
on  cyclic  plasticity-based  Fatigue  Indicator  Parameters  (FIPs)  that  reflect  probability  of  microstracture 
scale  crack  formation  and  microstructurally  small  crack  growth. 

For  most  applications,  large  numbers  of  experiments  are  necessary  to  meaningfully  quantify  any 
variability  in  fatigue  life  and  to  identify  any  change  in  the  mechanism  of  fatigue  damage  formation  as  a 
function  of  applied  loading  conditions  for  a  given  specimen  size.  Typically,  insufficient  experimental  data 
are  available  to  support  this  quantification.  Moreover,  the  mathematical  form  of  the  tails  of  the  probability 
distributions  for  the  driving  forces  for  fatigue  damage  formation,  such  as  the  local  distributions  of 
stress/strain,  are  not  well  characterized;  moreover,  the  dependence  of  the  character  of  those  tails  on  single 
and/or  interacting  microstracture  attributes  are  not  well  understood. 

In  this  work,  we  are  primarily  concerned  with  how  the  crystallographic  attributes  of  a  polycrystalline 
microstracture  (e.g.,  grain  orientation,  disorientation,  size,  and  shape  distributions)  affect  local  driving 
forces  for  fatigue  crack  formation  and  early  growth  in  HCF.  In  general,  plastic  strain  inhomogeneity  at  the 
grain  level  in  polycrystals  subjected  to  cyclic  loading  is  directly  linked  to  crystallographic  texture.  Winter 
et  al.  [7]  observed  that  plasticity  occurs  preferentially  in  grains  having  slip  systems  with  high  Schmid 
factors,  with  slip  localized  within  slip  bands.  Using  crystal  plasticity  simulations  to  calculate  distributions 
of  cyclic  slip  in  polycrystals,  Bennett  and  McDowell  [8]  demonstrated  that  distributions  of  slip  could  be 
quite  heterogeneous  in  HCF.  This  heterogeneity  is  directly  related  to  complex  interactions  between  grains 
of  differing  orientations.  Sauzay  and  Jourdan  [9]  explored  these  types  of  interactions  between  grains  by 
computationally  characterizing  the  distributions  of  elastic  stress  fields  around  small  grain  clusters  at  the 
free  surface  using  elastic  FE  simulations.  They  predicted  that  grain  interactions  could  affect  the  local 
resolved  shear  stress  by  as  much  as  18%  in  copper  and  austenitic  stainless  steels  depending  on  the  local 
orientations  of  the  neighboring  grains.  Inhomogeneity  of  elastic  stress  fields  corresponds  to  localization  of 
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plastic  strain  in  regions  of  stress  concentration  associated  with  the  jump  of  the  elastic  stiffness  across  grain 
boundaries  and  compatibility  requirements  of  the  polycrystal.  Specifically,  we  will  consider  the  effects  of 
local  crystallography  (e.g.,  phase,  grain  orientation,  grain  disorientation,  grain  topology,  etc.)  on  slip  in  a 
powder  metallurgy  (PM)  Ni-base  superalloy,  IN100. 

Ni-base  superalloys  are  predominantly  used  in  aircraft  gas  turbine  engines  due  to  their  high  strength 
and  creep  resistance  at  high  temperatures  that  is  conferred  by  coherent  y'  N'i  .AI  precipitates  of  Ll2  face- 
centered  cubic  (fee)  structure.  These  precipitates  are  dispersed  in  the  y  austenitic  Ni  solid  solution  matrix  of 
fee  crystal  structure  and  provide  excellent  resistance  to  slip.  Commonly,  fatigue  damage  formation  in 
polycrystalline  superalloys  has  been  linked  to  the  existence  of  large  pores  or  nonmetallic  inclusions 
introduced  during  processing.  Often,  inclusions  are  debonded  form  the  matrix  or  cracked  during  primary 
forming  processes.  During  loading,  the  stress  concentrations  at  inclusions/pores  often  lead  to  the  formation 
of  fatigue  cracks  [10-12].  However,  as  processing  techniques  improve,  cleaner  Ni-base  superalloys  are 
being  developed  that  have  lower  number  density  of  inclusions/pores;  consequently,  fatigue  cracks  are 
increasingly  observed  to  form  along  crystallographic  planes.  For  example,  Jha  et  al.  [13]  noted  that 
subsurface  fatigue  crack  formation  occurs  in  individual  grains  absent  of  any  voids/inclusions  in  Rene 
88DT,  particularly  at  lower  stress  amplitudes  ( i.e .,  in  the  HCF  regime).  In  this  case,  the  nucleation  region 
associated  with  the  size  of  the  crystallographic  facets  at  the  sites  of  fatigue  crack  formation  was  observed  to 
be  much  larger  than  the  average  grain  size.  This  suggests  that  cracks  tend  to  form  in  larger  grains.  In  the 
same  alloy,  Shyam  et  al.  [14]  noted  that  cracks  form  predominantly  in  larger  grains  or  at  inclusions  near 
large  grains.  Additionally  in  Rene  88DT,  Miao  et  al.  [15]  also  observed  that  most  critical  (i.e..  life  limiting) 
fatigue  cracks  in  the  HCF/VHCF  regime  initiate  crystallographically  away  from  the  surface.  The  grains  in 
which  these  cracks  form  were  observed  to  be  large  in  size  relative  to  the  average  grain  size.  Here  these 
grains  were  associated  with  higher  Schmid  factors,  indicating  they  are  oriented  favorably  for  slip.  Others 
have  shown  that  slip  bands  associated  with  shearing  of  the  second  phase  y'  precipitates  and  subsequent 
fatigue  damage  formation  have  been  observed  to  be  more  common  in  coarser  grained  superalloys,  while 
deformation  in  the  smaller  grains  has  been  observed  to  be  more  homogeneous  [16,  17]. 

The  analysis  of  the  slip  character  in  Ni-base  superalloys  has  provided  some  insight  into  fatigue  damage 
formation  in  these  materials  systems.  As  mentioned  previously,  the  y  austenitic  phase  in  these  Ni-base 


4 


Microstructure-Sensitive  Extreme  Value  Probabilities  for  High  Cycle  Fatigue  ofNi-Base  Superalloy  IN  100 

Craig  Przybyla,  David  L.  McDowell 

superalloys  has  a  fee  crystalline  lattice.  As  typical  in  fee  materials,  slip  is  expected  on  the  12  octahedral  {1 
1  1}<1  1  0>  slip  systems.  However,  contributions  of  an  additional  6  cube  {1  0  0}<1  1  0>  slip  systems  have 
been  observed  in  the  y  phase  in  many  Ni-base  superalloys  at  certain  elevated  temperatures  in  the  form  of 
macroscopic  slip  traces  [18].  Although  limited  cube  slip  is  possible  in  the  Ll2  ordered  y'  phase,  it  is  not 
commonly  considered  as  likely  to  occur  in  the  austenitic  y  phase.  Through  detailed  investigation  using 
transmission  electron  microscopy  (TEM),  Bettge  and  Osterle  [18],  have  postulated  that  this  cube  slip  along 
{10  0}  planes  is  actually  due  to  the  “zig-zag”  cross  slip  of  screw  dislocations  on  {1  1  1}  planes  when 
blocked  at  the  y/y’  interfaces,  as  depicted  in  Figure  1.  As  these  dislocation  migrate  within  the  channels 
between /'precipitates,  they  create  macroscopic  {10  0}  slip  traces. 

Cube  slip  has  been  observed  in  several  single  crystal  Ni-base  superalloys  over  a  wide  range  of 
temperatures.  For  example,  Miner  et  al.  [19]  investigated  slip  traces  in  single  crystal  of  Rene  N4  at  room 
temperature,  650°C,  760°C,  870°C,  and  980°C.  They  observed  slip  traces  along  {10  0}  planes  for  all  the 
temperatures  considered  when  the  crystals  were  oriented  with  the  {111}  slip  planes  perpendicular  to  the 
loading  direction.  Additionally,  at  increasingly  higher  temperatures,  {10  0}  slip  traces  were  observed  with 
increasing  frequency  in  the  crystals  oriented  in  directions  other  than  <1  1  1>  relative  to  the  loading 
direction.  Effectively,  cube  slip  accommodates  deformation  for  grains  in  “hard”  crystallographic 
orientations.  At  870°C  {111}  slip  traces  were  only  observed  for  crystals  oriented  near  [0  0  1]  and  [0  1  1] 
and  all  others  were  of  {1  0  0}  character.  Bettge  and  Osterle  [18]  observed  “zig-zag”  cube  slip  in  single 
crystals  of  SC  16  oriented  with  the  {111}  slip  planes  perpendicular  to  the  loading  direction.  Westbrooke  et 
al.  [20]  observed  slip  bands  traces  along  {100}  planes  in  single  crystals  of  a  Pratt  &  Whitney  alloy  for  a 
crystal  with  unfavorable  orientation  relative  to  slip  on  {1  1  1}  planes  and  octahedral  slip  for  a  crystal 
oriented  for  loading  in  the  |”l  10J  direction.  They  also  observed  slip  traces  for  the  crystal  tested  in  the  [1  0 

0]  direction  that  were  neither  of  [1  1  1}  or  {1  0  0}  character.  The  morphology  and  distribution  of  y  ’ 
appears  to  dictate  the  character  of  slip  in  these  types  of  superalloys.  For  example,  Nitz  et  al.  [21]  examined 
slip  in  single  crystals  of  NIMONIC  105  at  temperatures  from  10°C  to  877°C  via  TEM  analysis.  In  this 
material  they  only  observed  octahedral  slip.  We  note  that  in  this  particular  material  the  y’  precipitates  are 
small  relative  to  the  other  materials  discussed  here  (~17  nm)  and  are  spherical  in  their  morphology  whereas 
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in  most  single  crystals  the  y’  precipitates  are  an  order  of  magnitude  larger  and  tend  to  have  a  cuboidal 
morphology.  The  slip  character  of  these  single  crystals  is  summarized  in  Table  1. 

Various  slip  modes  have  also  been  observed  in  fatigue  tests  of  various  polycrystalline  Ni-base 
superalloys.  The  slip  behavior  of  various  polycrystalline  alloys  is  summarized  in  Table  2.  Generally  in 
these  alloys  fatigue  cracking  propagates  along  {111}  planes  at  lower  homologous  temperatures  and  along 
{10  0}  planes  at  moderate  temperatures  (i.e.,  between  500°C  to  800°C).  We  point  out  that  this  is  what  was 
observed  in  fatigue  testing  in  IN  100  [22].  In  contrast,  when  the  y’  precipitates  are  more  spherical  and 
smaller  in  size  and  have  lower  volume  fraction,  slip  and  cracking  along  {1  0  0}  planes  has  not  been 
observed,  such  as  in  Rene  88  DT  [23]. 

Statistical  treatments  of  fatigue  have  primarily  been  based  on  large  numbers  of  experiments.  Typically, 
one  observes  a  wide  range  of  scatter  for  a  given  component  in  the  number  of  cycles  to  failure,  particularly 
in  HCF.  For  example.  In  Figure  2  we  can  see  significant  variability  in  the  overall  fatigue  lives  in  INI 00 
across  a  range  of  applied  stress  magnitudes  [24].  The  variability  of  fatigue  life  is  assessed  by  extensive 
experimentation  to  obtain  a  statistically  significant  sample.  Designers  then  use  these  data  to  predict 
component  life  with  an  acceptable  level  of  risk.  Such  data  collection  requires  significant  time  and  resources 
and  does  not  necessarily  provide  understanding  of  the  mechanism(s)  involved  in  dictating  variability.  In 
addition,  the  resulting  predictions  often  change  with  sample  size.  The  problem  is  further  complicated  when 
multiple  mechanisms  of  damage  formation  are  observed  [13,  24,  25].  In  these  cases,  multiple  competing 
mechanisms  of  fatigue  damage  formation  and  fatigue  crack  growth  make  it  difficult  to  accurately  predict 
fatigue  life  because  different  mechanisms  operate  at  different  stress  magnitudes  and  often  microstructure 
attributes  vary  spatially  in  the  material.  For  example,  multiple  modes  of  failure  are  evident  in  the 
cumulative  probability  plot  in  Figure  2(b)  at  the  stress  magnitude  of  1150MPa.  In  this  case,  the  authors 
assert  that  the  bi-modal  nature  of  the  failure  probability  relates  to  the  fact  that  there  exist  two  different 
competing  failure  mechanisms  [24]. 

Processes  of  HCF  crack  formation  and  early  growth  depend  on  a  few  key  attributes  (or  sets  of 
attributes)  such  as  the  largest  inclusion  or  void  within  an  entire  population  of  such  features.  There  is 
significant  scatter  in  experimentally  measured  HCF  life  because  attributes  that  govern  failure  lie  within  the 
tails  of  joint  distributions  of  attributes  and  responses,  necessitating  consideration  of  extreme  value  statistics. 
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For  example,  if  the  key  feature  of  the  fatigue  life  in  a  particular  material  is  determined  by  the  largest 
inclusion  of  a  certain  type,  an  understanding  of  the  distribution  of  inclusion  sizes  in  the  critically  stressed 
regions  would  be  essential  to  modeling  the  F1CF  life  for  that  particular  material. 

The  problem  of  the  largest  inclusion  size  in  clean  steels  has  been  considered  by  several  different 
workers  and  is  summarized  in  a  review  by  Atkinson  [26].  In  this  case,  different  models  exist  that  predict 
the  fatigue  resistance  of  these  steels  based  on  the  largest  inclusion  size.  Specifically,  these  models  predict 
the  maximum  inclusion  size  using  an  assumed  log-normal  distribution,  statistics  based  on  the  extreme  value 
distribution  function  of  Gumbel  [27],  or  a  statistical  method  based  on  the  Generalized  Pareto  Distribution 
(GPD).  All  of  these  methods,  however,  assume  failure  based  on  a  single  attribute  (e.g.,  inclusion  size)  and 
do  not  consider  how  interacting  attributes  affect  the  damage  processes  of  interest.  For  example,  in  many 
cases  the  driving  forces  for  damage  formation  around  an  inclusion  might  depend  as  much  on  the  orientation 
of  the  grains  that  affect  the  local  slip  processes  in  the  matrix  around  an  inclusion  as  on  the  size  of  the 
inclusion  itself.  We  also  argue  that  a  framework  that  only  considers  purely  geometric  attributes  such  as 
inclusion  size,  grain  size,  or  grain  orientation  is  in  general  insufficient  to  quantify  the  scatter  in  F1CF  life. 
The  additional  attribute(s)  of  material  response  (e.g.,  stress,  strain,  plastic  strain  and  other  driving  force 
parameters)  coupled  with  geometric  attributes  of  microstructure  facilitate  pursuit  of  joint  statistics  of 
extreme  value  type  that  are  relevant  to  minimum  fatigue  life  design. 

In  the  present  work,  we  propose  a  new  statistical  framework  that  links  extreme  value  probabilities  of 
various  response  parameters  to  the  correlated  microstructure  attributes  associated  with  the  extreme  value 
response.  Here  we  exercise  this  statistical  framework  using  computational  simulations  that  estimate  the 
fatigue  response  in  a  Ni-base  superalloy  IN100  using  a  crystal  plasticity  formulation  via  the  finite  element 
(FEM)  method.  Finally,  we  compare  our  results  to  a  limited  number  of  fatigue  experiments  available  in  the 
literature  for  this  particular  material  system. 

2.  Methodology 

2.1  Extreme  Value  Marked  Radial  Microstructure  Correlation  Functions 

The  objective  is  to  develop  a  statistical  framework  that  is  able  to  quantify  the  coupling  of  the  extreme 

values  of  certain  FIPs  to  the  key  microstructure  attributes  that  are  associated  with  these  extreme  value 


7 


Microstructure-Sensitive  Extreme  Value  Probabilities  for  High  Cycle  Fatigue  ofNi-Base  Superalloy  IN  100 

Craig  Przybyla,  David  L.  McDowell 


response  parameters.  To  achieve  this  objective  we  postulate  that  our  proposed  statistical  framework  must 
meet  three  primary  requirements: 

1.  It  must  contain  information  regarding  both  the  distributions  of  microstructure  attributes  and 
potential  interactions  or  correlations  between  multiple  microstructure  attributes. 

2.  It  must  link  the  distributions  of  microstructure  attributes  and  correlations  between  various 
microstructure  attributes  to  distributions  of  local  response  parameters  (e.g.,  FIPs). 

3.  It  must  be  able  to  address  the  extreme  value  nature  of  fatigue  crack  formation  and  early  growth 
scenarios. 

We  propose  a  two  component  statistical  framework  to  satisfy  these  requirements.  Given  a  spatial 
window  Q,  we  define  the  probability  distribution  of  the  extreme  value  response  parameter  a  as  Fex  (a|Q)> 

which  is  associated  with  the  probability  that  the  response  parameter  of  value  a  is  the  extreme  value  for  a 
sampled  statistical  volume  element  of  microstructure  O.  This  follows  the  ideas  of  extreme  value  (or 
importance)  sampling  from  classical  extreme  value  statistics  of  Gumbel  [27,  28].  Corresponding  to  the 
probability  distribution  of  the  extreme  value  response  parameter  is  the  extreme  value  marked  radial 
correlation  function,  ^maxM  \  r,Cl) ,  which  is  associated  with  the  probability  of  finding  of  a  sphere 

centered  at  the  microstructure  attribute  ji  corresponding  to  the  location  of  the  maximum  response  parameter 
a  in  the  microstructure  window  O,  with  microstructure  attribute  /?'  at  a  distance  within  r  to  r  +  dr  from  /?  in 
any  direction.  Although  based  on  the  notion  of  the  marked  correlation  functions  introduced  by  Pryz  [29], 
this  new  construct  is  fundamentally  unique  in  its  framing  of  extreme  value  statistics.  Here  we  have  based 
the  extreme  value  marked  correlation  function  on  the  classical  radial  distribution  function,  but  it  can  just  as 
easily  be  constructed  using  n-point  correlation  functions,  lineal  path  correlation  functions,  etc. 

This  two  component  statistical  construct  captures  both  the  extreme  value  response  of  the 
microstructure  as  represented  by  the  response  parameter  a  and  the  biased  correlations  of  microstructure 
attributes  between  ji  and  ji'  in  the  neighborhood  of  the  observed  extreme  values  of  a  in  a  microstructure 
window  O.  In  this  way,  spatial  correlations  are  identified  between  microstructure  attributes  that  have  a  high 
probability  of  existing  in  the  neighborhood  of  an  extreme  value  response  parameter.  This  sampling  is 
performed  on  microstructure  window  Q  and  is  expected  to  depend  on  the  size  of  Q.  The  method  requires 
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multiple  instantiations  of  Q  to  effectively  sample  both  F“(a|Q)  and  i?max(a)  (/?,/?' |  r,Q)  to  build  up  the 

tail  of  the  probability  distribution  corresponding  to  extreme  value  response  neighborhoods.  The  present 
framework  facilitates  comparison  among  multiple  microstructures  for  a  given  material  or  comparison  of 
significant  microstructure  attributed  being  considered  in  the  sampling  different  materials  which  can  support 
materials  design. 

Constructing  the  extreme  value  distribution  with  the  corresponding  extreme  value  marked  correlation 
function  requires  a  significant  number  of  simulations/experiments  if  they  are  to  be  considered  statistically 
meaningful.  This  can  require  extensive  processing  time  for  both  simulations  and  data  analysis  depending  on 
the  complexity  of  the  models  being  analyzed.  The  number  of  samples  required  depends  on  the  response 
parameter  considered  and  material  being  analyzed.  Analogous  experimental  data  are  not  typically  available 
and  are  very  expensive  to  obtain.  Typical  fatigue  experiments  will  not  provide  enough  information  to 
construct  these  coupled  distribution  functions  because  such  experiments  consist  of  only  a  few  data  points  at 
each  stress/strain  amplitude  considered.  In  addition,  experimentally,  fatigue  hot  spots  are  identified  after 
failure  and  characterization  of  the  sites  of  fatigue  damage  formation  and  their  neighborhoods  is  very 
cumbersome  and  typically  requires  destructive  sectioning.  Thus,  although  this  proposed  framework  is 
general  in  its  application  to  experiments  and/or  simulations,  we  envision  that  it  will  be  much  more 
practically  applied  to  simulations  that  are  substantiated  with  limited  experimentation. 

2.2.  Simulation  Strategy 

To  exercise  the  foregoing  extreme  value  statistical  framework  for  polycrystalline  Ni-base  superalloys, 
we  choose  a  strategy  that  estimates  the  proposed  statistical  parameters  via  the  simulation  of  multiple 
microstructure  instantiations.  Specifically,  we  examine  the  influence  of  polycrystalline  microstructure  on 
probability  of  crystallographic  fatigue  crack  formation  in  the  Ni-base  superalloy  IN100.  We  recognize  that 
fatigue  damage  formation  in  this  particular  material  system  primarily  relates  to  the  presence  of  non-metallic 
inclusions  and/or  pores.  However,  we  have  noted  that  Ni-base  superalloys  have  been  developed  with  much 
smaller  populations  of  these  types  of  inclusions/voids  that  do  fail  predominantly  crystallographically  (e.g., 
[23]).  Thus,  in  this  work  we  are  exploring  how  this  variant  of  IN100  would  fail  crystallographically  in  the 
absence  of  non-metallic  inclusions  or  voids.  In  addition,  the  improved  understanding  of  extreme  value 
behavior  of  the  orientation  dependent  slip  processes  in  IN  100  developed  here  can  later  be  coupled  with  the 
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existence  of  inclusions/voids  to  better  current  predictive  models  of  fatigue  crack  formation  in  this  material 
system.  To  accomplish  this  objective,  digital  representations  of  microstructures  are  simulated  via  the  FE 
package  ABAQUS  [30]  coupled  with  a  microstructure-sensitive  crystal  plasticity  model  for  IN  100  [31,  32] 

We  briefly  describe  the  crystal  plasticity  model  for  IN  100.  Generally  superalloys  are  more  complex  to 
model  because  of  their  tension-compression  asymmetry  and  the  non-Schmid  characteristics  of  the  Ni3Al 
phase  due  to  dislocation  core  spreading  effects.  The  grain  size  of  the  primary  y  phase  along  with  the  size 
and  spacing  of  the  coherent  y’  precipitates  can  greatly  affect  the  material  response  of  these  material 
systems.  For  example,  dislocation  mechanisms  change  from  precipitate  shearing  to  bypassing  by  Orowan 
looping  depending  on  the  spacing,  size  and  volume  fraction  of  the  y  ’  precipitates,  as  well  as  the  amplitude 
of  applied  strain.  In  this  model  [31,  32],  the  precipitates  are  not  explicitly  modeled,  but  the  influence  of 
primary,  secondary  and  tertiary  y  ’  precipitates  are  included  through  certain  internal  state  variables.  Also 
incorporated  in  the  model  are  the  average  effect  of  the  grain  size  on  the  constitutive  response  and  the 
dependence  of  the  flow  stress  on  crystallographic  orientation.  Shenoy  et  al.  [31,  32]  consider  both  the 
standard  12  octahedral  <1  1  0>{1  11}  slip  systems  and  6  cube  slip  systems  {10  0}<1  1  0>.  As  previously 
mentioned,  the  cube  slip  systems  describe  “zig-zag”  octahedral  slip  that  has  been  observed  at  the  y-y' 
interface  in  the  y  matrix  and  cube  slip  that  occurs  in  y'  precipitates  at  higher  homologous  temperatures  for 
crystals  in  “hard”  orientations  for  octahedral  slip.  This  particular  model  has  been  calibrated  for  IN100  at 
650°C. 

2.3  Digital  Microstructure  Models 

Using  a  crystal  plasticity  formulation  allows  these  digital  microstructure  volumes  to  be  statistical 
volume  elements  (SVEs)  for  local  plasticity  response  (a  key  driving  force  for  fatigue  damage  formation). 
Each  statistical  volume  element  for  a  desired  response  (i.e.,  grain  scale  plasticity)  is  designed  to  be  large 
enough  such  that  the  response  parameter  of  interest  at  a  particular  location  is  unaffected  by  statistical 
variations  in  the  microstructure  at  distances  on  the  order  of  the  size  of  the  volume  (i.e.,  correlation  length  is 
less  than  size  of  SVE);  however,  the  volume  is  not  so  large  that  it  contains  a  statistically  representative  set 
of  responses  for  that  particular  response  parameter.  This  is  in  contrast  to  a  representative  volume  element 
(RVE)  that  is  defined  such  that  the  distribution  of  local  response  or  effective  response  of  the  volume  will 
not  change  based  on  where  the  RVE  is  sampled  from  the  microstructure  ensemble  (i.e.,  bulk  material),  or  if 
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it  is  further  increased  in  size.  In  some  cases,  information  obtained  from  multiple  SVEs  can  constitute  a 
RVE  level  characterization.  We  note,  however,  that  it  is  sometimes  the  case  that  for  a  particular  response  a 
RVE  of  reasonable  size  is  unattainable.  This  is  particularly  true  when  trying  to  capture  extreme  value 
distributions  of  microstructure  response  that  affect  rare  event  damage  formation  and  growth.  This  implies 
that  complete  characterization  of  the  tail  of  the  probability  distribution  function  (PDF)  for  fatigue  response 
requires  a  very  large  RVE. 

In  the  present  case,  we  seek  to  simulate  multiple  SVEs  of  sufficient  size  such  that  the  lower  order 
moments  of  the  cyclic  plasticity  response  on  the  scale  of  the  grains  are  unaffected  by  further  increasing  the 
volume  of  the  SVE.  Of  course,  the  higher  order  moments  affect  the  tail  of  the  PDF  and  the  RVE  size  for 
invariance  of  these  moments  would  be  quite  large,  in  general.  We  also  desire  that  the  SVE  be  suitably  large 
to  serve  as  a  RVE  in  terms  of  the  overall  effective  (i.e.,  macroscopic)  elastic/plastic  response.  In  other 
words,  we  want  the  stress-strain  response  taken  from  any  microstructure  instantiation  to  be  representative 
of  any  other.  The  RVE  for  this  effective  response  is  insensitive  to  the  higher  order  moments  of  the  slip 
distribution,  and  is  not  the  same  as  a  RVE  for  the  extreme  value  distribution  of  FIPs.  This  construct  follows 
the  same  ideas  behind  SERVEs  as  defined  by  Swaminathan  and  Ghosh  [33,  34].  Being  able  to  simulate 
variations  in  the  local  response  via  these  types  of  SVEs  enables  us  to  computationally  explore  how 
microstructure  heterogeneity  affects  fatigue  variability. 

A  Voronoi  tessellation  is  a  mathematical  construct  that  describes  a  space  partitioned  around  predefined 
centers  such  that  all  points  closer  to  a  given  center  than  any  other  is  associated  with  that  center.  In  this 
manner,  the  space  is  divided  up  into  a  set  of  space  filling  convex  polyhedral.  Voronoi  tessellations  have 
been  used  in  many  different  applications  and  the  geometrical  concept  behind  this  construction  is  well 
defined  in  the  literature  [35].  Fast  Voronoi  generation  algorithms  are  readily  available  [36].  Goss  and  Li 
[37]  evaluated  Voronoi  tessellations  based  on  structure,  topology  and  statistics  relative  to  typical 
polycrystalline  materials  and  observed  that  both  topology  and  various  statistical  properties  do  not  agree 
well  with  those  measured  experimentally.  They  suggest  that  more  realistic  Voronoi  tessellated  models  for 
polycrystalline  materials  can  be  generated  using  optimization  methods  to  randomly  perturb  the  centroids  of 
the  Voronoi  tessellated  cells  to  match  certain  statistical  distributions  measured  from  actual  material  systems 
like  grain  size,  grain  volume,  grain  boundary  length,  etc.  It  was  also  suggested  that  optimization  methods 


11 


Microstructure-Sensitive  Extreme  Value  Probabilities  for  High  Cycle  Fatigue  ofNi-Base  Superalloy  IN  100 

Craig  Przybyla,  David  L.  McDowell 

could  be  used  to  match  various  texture  measurements  such  as  orientation  or  misorientation/disorientation 
(cf.  [38]).  Modifying  random  Voronoi  tessellations  to  match  more  realistic  distributions  of  grain  size  has 
also  been  performed  earlier  by  others  [39,  40]  to  create  RVEs  for  FE  simulations.  Recently,  Zhang  et  al. 
[41,  42]  used  modified  Voronoi  tessellations  that  were  optimized  with  a  simulated  annealing  algorithm  to 
fit  experimentally  characterized  statistics  of  orientation,  disorientation,  and  phase  volume  fraction  in 
multiphase  Ti-6A1-4V  to  construct  FE  models  for  crystal  plasticity  simulations.  Additionally,  Shenoy  et  al. 
[31]  constructed  digital  microstructure  instantiations  for  IN  100  using  the  same  methods.  Similar  algorithms 
are  used  here  to  fit  distributions  of  grain  size,  orientation  and  disorientation  to  construct  SVEs  for  IN  100.  A 
SVE  generated  for  IN  100  along  with  an  imposed  voxellated  mesh  can  be  seen  in  Figure  3.  The  target  and 
optimized  distributions  for  the  grain  volume  normalized  by  the  target  average  grain  volume  are  given  in 
Figure  4. 

The  microstructures  simulated  in  this  work  are  constructed  via  a  C++  computer  program  that  generates 
these  statistical  volume  elements  of  realistic  polycrystalline  microstructures.  The  actual  Voronoi 
tessellation  is  computed  using  QFIULL  [36]  (see  www.qhull.org).  The  simulated  polycrystalline 
microstmctures  are  generated  using  this  modified  Voronoi  tessellation  algorithm  such  that  all  points  closer 
a  grain  center  than  any  other  grain  center,  belong  to  that  grain.  An  additional  parameter  is  also  input  by  the 
user  that  defines  the  ratio  of  the  minimum  distance  between  any  two  grain  centers  to  average  distance 
between  the  grain  centers.  This  prevents  grain  centers  lying  too  close,  thus  also  preventing  the  generation 
of  very  small  voronoi  cells  (i.e.,  grains).  Fitting  the  simulated  distributions  of  grain  size,  grain  orientation, 
and  grain  disorientation  is  performed  in  the  same  manner  as  demonstrated  previously  by  others  [31,41,  42]. 
In  our  case,  we  chose  a  random  orientation  distribution  function  (as  commonly  observed  in  IN  100)  and  a 
random  Mackenzie  disorientation  distribution  function  [43,  44]. 

Once  a  microstructure  volume  or  instantiation  is  constructed  and  optimized  to  fit  the  target  distributions 
of  grains  size,  grain  orientation,  and  grain  disorientation,  the  program  writes  the  ABAQUS  input  files  that 
contain  all  the  necessary  information  to  simulate  the  microstructures.  Information  in  these  input  files 
includes  the  details  of  the  geometry,  mesh  and  boundary  and  loading  conditions.  The  initial  simulation  cell 
geometry  is  defined  by  the  user  defined  values  of  edge  length  and  mesh  density  (the  number  of  elements 
along  each  edge  of  the  simulated  block  of  microstructure).  The  user  also  defines  the  target  distributions  for 
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the  grain  size,  grain  orientation  and  grain  disorientation.  Once  the  global  simulation  cell  geometry  is 
generated  (with  grain  centers  defined),  the  elements  and  nodes  are  created  based  on  the  mesh  density  input 
by  the  user.  Using  the  known  edge  length  and  mesh  density,  the  element  edge  length  and  total  size  is 
calculated.  The  elements  are  constructed  as  an  array  of  cuboidal  elements  (type  C3D8R  in  ABAQUS)  each 
with  identical  dimensions.  These  particular  elements  employ  reduced  integration  to  speed  up  computation 
time.  Cuboidal  elements  are  also  convenient  because  the  overall  simulation  cell  was  assumed  to  be  cubic  in 
this  model.  In  addition,  a  voxellated  mesh  greatly  simplifies  the  application  of  periodic  boundary 
conditions  in  all  directions. 


2.4  Microstructure  Attributes  and  Fatigue  Indicator  Parameters 

The  strain-life  approach  in  fatigue  has  been  applied  mainly  in  the  low  cycle  and  transition  fatigue 

regimes;  however,  microplasticity  (i.e.,  heterogeneous  plasticity  at  the  scale  of  microstructure)  is  relevant 

to  understanding  fatigue  crack  formation  and  early  crack  propagation  in  HCF  [2,  4],  Coffin  [45]  and 

Manson  [46]  independently  proposed  a  power  law  correlation  of  plastic  strain  amplitude,  As:.p/2,  with  the 

number  of  cycles,  Nf,  to  fatigue  crack  initiation  (formation  plus  small  crack  growth  to  some  predefined 

size)  for  laboratory  specimens,  i.e., 

where  s'  is  the  fatigue  ductility  coefficient  and  c  is  the  fatigue  ductility  exponent  [47]. 

Two-parameter  approaches  have  been  developed  as  a  multiaxial  generalization  of  plastic  strain  amplitude 
in  an  extension  of  Eq.  (1).  Critical  plane  plastic  strain  range-based  approaches  have  been  introduced  that 
consider  cyclic  plastic  shear  strain  as  the  primary  driving  force  for  fatigue  crack  formation  and  early 
growth,  with  a  modifying  influence  of  normal  stress  or  strain  acting  on  the  plane  of  maximum  reversed 
plastic  shear  strain.  In  the  present  work  we  regard  such  quantities  as  Fatigue  Indicator  Parameters  (FIPs)  in 
view  of  their  relation  to  probability  of  crack  formation  and  early  growth.  Fatemi  and  Socie  [48]  and  Fatemi 
and  Kurath  [49]  proposed  a  parameter  for  shear-dominated  crack  initiation  that  accounts  for  the  observation 
that  fatigue  cracks  initiate  on  planes  of  maximum  shear  for  multiaxial  loading  conditions,  i.e., 
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where  A;/(’,ix  is  the  maximum  range  of  plastic  shear  strain  and  cr"  is  the  maximum  normal  stress  acting  on 

this  plane.  The  stress,  cr"ax ,  accounts  for  the  effect  of  normal  stress  to  the  plane  of  small  crack  formation 

and  growth,  and  is  weighted  by  the  material  constant  k.  This  term  is  normalized  by  the  yield  strength. 
Larger  magnitudes  of  PFS  indicate  a  higher  driving  force  for  fatigue  damage.  The  Fatemi-Socie  (FS) 
parameter  has  been  used  to  correlate  fatigue  damage  formation  over  a  large  number  of  grains  for  a  range  of 
multiaxial  loading  conditions  for  materials  with  extended  Stage  I  dominant  regions  such  as  1 045  steel  and 
IN718  [50,  51].  Although  the  Fatemi-Socie  FIP  was  originally  only  applied  at  the  scale  of  laboratory 
specimens,  several  workers  have  applied  this  and  similar  FIPs  at  the  microscale  using  crystal  plasticity 
simulations.  For  example,  Bennett  et  al.  [52]  correlated  microslip  and  mixed-mode  behavior  of 
microstmrally  small  crack  growth  to  the  Fatemi-Socie  FIP  and  considered  its  distribution  over  the  grains; 
moreover,  the  Fatemi-Socie  FIP  distributes  heterogeneously  over  the  grains.  Moreover,  they  also  observed 
that  the  Fatemi-Socie  FIP  correlates  well  with  distributions  of  small  fatigue  cracks  within  the  grains. 
Bennett  et  al.  [8]  also  used  the  Fatemi-Socie  FIP  to  explore  local  grain  orientation  distribution  effects  on 
microslip  in  HCF.  Other  workers  [53,  54]  observed  good  agreement  between  the  locally  applied  Fatemi- 
Socie  FIP  (determined  using  crystal  plasticity  calculations)  and  microstructurally  small  crack  growth  within 
the  first  few  grains  based  on  experiments  in  structural  steel  S460N.  Additionally,  Dunne  et  al.  [55]  used 
detailed  simulations  coupled  with  experiments  to  show  that  crystal  plasticity  can  correlate  precisely  to 
localized  plastic  slip  (and  subsequent  fatigue  crack  formation)  in  realistic  Ni-base  superalloy  and  Ti  alloy 
microstructures.  More  recently  Findley  and  Saxena  [56]  used  the  Fatemi-Socie  FIP  to  examine  the  effect 
microstructure  attributes  such  as  grain  size  have  on  the  local  driving  forces  for  fatigue  damage  formation. 
Such  work  demonstrates  the  ability  of  locally  applied  FIPs  determined  via  crystal  plasticity  simulations  to 
correlate  microstructure  scale  slip  with  fatigue  crack  formation  and  early  stages  of  microstructurally  small 
crack  growth.  In  essence,  there  is  similitude  in  scaling  the  relations  between  Pls  and  the  number  of  cycles 

to  form  and  growth  a  fatigue  crack  to  different  length  scales  within  the  range  of  several  microns  to  several 
hundred  microns. 

As  discussed  by  McDowell  [2],  these  types  of  critical  plane  approaches  for  shear-dominated 
microcracking  are  related  to  Stage  I  shear-dominated  propagation  of  MSCs  with  an  influence  of  the  normal 
stress  to  the  slip/crack  plane.  Early  work  by  Hoshide  and  Socie  [57]  considered  the  affect  of  microplastic 
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strain  on  fatigue  damage  formation  and  MSC  crack  growth  under  mixed-mode  (I-II)  loading.  They 
correlated  A J  from  EPFM  via  a  Paris  type  law  with  fatigue  crack  growth  behavior  of  small  cracks. 
Moreover,  Hoshide  and  Socie  [57]  related  A J  to  the  ranges  of  crack  opening  displacement  (COD)  and/or 
crack  sliding  displacement  (CSD)  for  mode  I  and  mode  II,  respectively.  Later,  McDowell  and  Berard  [50] 
employed  A J  and  argued  for  similarity  of  A CTD  with  the  Fatemi-Socie  parameter,  which  assists  in 
understanding  the  aforementioned  notion  of  similitude. 

3.  Results 

3.1  Mesh  Quality  Study 

As  seen  in  Figure  3,  we  impose  a  voxellated  mesh  to  model  the  generated  microstructures.  A  detailed 
study  was  performed  to  investigate  the  sensitivity  of  the  convergence  of  the  Fatemi-Socie  FIP  to  these 
types  of  voxellated  meshes.  The  parameters  for  the  microstructures  that  were  generated  using  the  Voronoi 
based  microstructure  generator  are  given  in  Table  3.  Four  different  microstructures  were  generated  using 
the  Voronoi  based  microstructure  generator  for  this  mesh  study  and  are  hereafter  referred  to  as 
microstructures  A,  B,  C  and  D.  Two  of  the  meshes  with  different  numbers  of  elements  for  Microstructure  C 
are  shown  in  Figure  5  with  contours  of  accumulated  plastic  strain  after  three  cycles.  In  each  case,  the 
elements  belonging  to  the  grain  with  the  extreme  value  grain  averaged  Fatemi-Socie  parameter  of  highest 
magnitude  among  all  of  the  grains  in  the  volume  is  highlighted  in  red.  In  the  case  of  Microstructure  C  as 
seen  in  Figure  5,  the  grain  highlighted  in  red  is  contained  entirely  in  the  interior  of  the  volume  (i.e.  the 
grain  does  not  intersect  the  boundary)  except  for  one  element  on  the  x-y  surface.  A  description  of  the 
different  meshes  considered  is  given  in  Table  4.  For  the  simulations,  periodicity  was  imposed  in  all 
directions.  The  simulations  were  cycled  in  strain  control  with  a  maximum  strain  of  0.5%  under  completely 
reversed  loading  (i.e.,  R=-l).  The  distribution  of  the  grain  volume  normalized  by  the  target  average  grain 
volume  of  1.0  x  10"6mm3  was  fit  to  a  log-normal  distribution  with  the  mean  and  standard  deviation  of  the 
natural  logarithm  of  the  normalized  grain  volume  given  by  -14.0  and  0.7,  respectively.  A  strain  rate  of 
0.002s"1  was  employed  in  these  simulations.  All  calculations  were  performed  over  the  entire  volume. 

The  convergence  of  the  extreme  value  Fatemi-Socie  FIP  as  calculated  over  a  single  element  and  over  all 
the  elements  of  each  grain  is  given  in  Figure  6.  In  Figure  6(a),  we  can  see  that  the  magnitude  of  the  extreme 
value  Fatemi-Socie  FIP  changes  by  as  much  at  20%  across  the  meshes  with  20  to  28  elements  along  an 
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edge  and  convergence  is  not  achieved.  However,  with  grain  size  averaging  volumes,  variation  in  the 
extreme  value  FIP  reduces  to  less  than  5%  in  most  cases  across  meshes  with  14  to  28  elements  along  an 
edge.  The  grain  with  the  observed  extreme  value  FIP  did  not  change  when  there  were  12  or  more  elements 
along  an  edge  or  when  the  ratio  of  the  element  volume  to  average  grain  volume  was  greater  than  or  equal  to 
0.064. 

3.2  Statistical  Volume  Elements  and  Size  for  IN100 

As  mentioned  previously,  we  define  our  SVEs  such  that  the  local  plastic  response  on  the  scale  of  the 

grains  is  unaffected  by  further  increasing  the  volume  of  the  SVE.  Additionally,  we  require  that  the  overall 

effective  (i.e.,  macroscopic)  elastic/plastic  response  is  nearly  the  same  for  each  SVE.  To  quantify  to 

influence  of  SVE  size  on  the  FIPs  we  simulated  several  SVEs  of  different  size.  Specifically,  we  constructed 

several  instantiations  of  cuboidal  arrays  of  randomly  oriented  grains.  These  cuboidal  grain  arrays  consist  of 

blocks  with  3,  5,  7,  and  9  grains,  respectively,  along  each  edge  with  random  periodic  boundary  conditions 

applied  in  all  directions.  In  this  manner  the  center  grain  has  1,  2,  3,  and  4  neighbors,  respectively,  in  each 

direction  between  itself  and  the  applied  periodic  boundary  conditions.  In  these  simulations,  these  SVEs 

were  cycled  one  time  under  completely  reversed  loading  (i.e.,  R=-l)  with  a  maximum  strain  magnitude  of 

1.0%.  The  quasistatic  strain  rate  of  0.004s'1  was  used  in  these  simulations.  The  grains  were  dimensioned  9 

pm  along  each  edge  and  consist  of  a  total  of  27  quadrilateral  elements  with  reduced  integration  (type 

C3D8R  in  ABAQUS)  or  3  elements  along  each  edge.  The  orientation  of  the  individual  grains  was  random 

except  for  the  center  grain  which  was  oriented  in  the  <0  0  1>,  <2  1  25>,  and  <1  1  1>  directions, 

respectively.  We  can  see  in  Figure  8  that  crystals  oriented  near  the  <2  1  25>  and  <1  1  1>  directions  exhibit 

the  maximum  apparent  Schmid  factors  for  octahedral  and  cuboidal  slip,  respectively.  The  Fatemi-Socie  FIP 

averaged  over  the  differently  oriented  center  grains  for  the  differently  sized  SVEs  are  given  in  Figure  7.  In 

all  cases,  there  is  a  significant  change  in  the  grain  averaged  FIP  of  the  center  grain  when  the  SVE  is 

increased  from  3  to  5  grains  along  each  edge.  However,  as  the  number  of  grains  along  each  edge  are 

increased  further  from  4  to  5,  5  to  7,  or  7  to  9,  further  changes  grain  averaged  FIP  of  the  center  grain  are 

much  less.  Thus,  for  the  purposes  of  this  work,  we  assume  that  the  primary  zone  of  influence  includes  the 

two  nearest  neighbors  of  any  one  grain  and  that  grains  further  out  than  that  have  much  less  influence.  The 

SVEs  for  our  microstructures  for  the  subsequent  study  are  scaled  accordingly. 
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3.3  Extreme  Value  Fatigue  Response  and  Microstructure  in  IN100 

To  investigate  the  influence  of  grain  orientation,  disorientation  and  grain  size  on  IN  100,  multiple  SVEs 

have  been  constructed  using  the  Voronoi  tessellation  based  microstructure  generator  and  simulated  via 
FEM,  as  previously  described.  Multiple  SVEs  for  the  local  plastic  response  were  generated  and  strained 
uniaxially  at  0.3%,  0.5%  and  0.7%  strain  at  a  quasistatic  strain  rate  of  -0.002s"1.  The  macroscopic  yield  is 
observed  for  this  particular  IN  100  at  approximately  near  1.0%  strain.  The  cyclic  strains  were  completely 
reversed  (i.e.,  R  =  -1).  Each  simulation  was  cycled  three  times  to  allow  for  an  initial  shakedown  period 
before  the  results  were  analyzed.  Periodic  boundary  conditions  were  applied  in  all  directions.  As  in  the 
previously  described  mesh  study  and  SVE  size  study,  the  voxellated  meshes  of  these  generated 
microstructures  consist  of  quadrilateral  elements  with  reduced  integration  (type  C3D8R  in  ABAQUS).  The 
number  of  grains  in  the  simulations  was  allowed  to  vary  to  best  fit  the  imposed  grain  size  distribution  and 
ranged  between  325-375  grains.  Each  SVE  was  dimensioned  to  0.150  mm  along  each  edge.  The 
distribution  of  grain  volumes  normalized  by  the  target  average  grain  volume  of  8.0xl0"6mm3  was  fit  to  a 
log  normal  distribution  with  a  mean  of  -12  and  standard  deviation  of  0.4.  This  corresponds  to  a  target 
average  cube  root  grain  size  of  -0. 020mm.  There  are  24  elements  along  each  edge  of  the  SVE  or  13824 
elements  in  all.  The  ratio  of  the  element  volume  to  the  average  target  grain  volume  is  0.030.  We  simulated 
a  total  of  50  SVEs  cycled  at  0.3%  strain,  100  SVEs  cycled  at  0.5%  and  50  SVEs  cycled  0.7%  strain. 

To  evaluate  the  convergence  of  the  extreme  value  distribution  {i.e.,  Fex  (a  |  Q) )  of  the  Fatemi-Socie  FIP, 

we  considered  the  shape  of  the  distribution  as  determined  from  sample  sizes  of  25,  50,  75  and  100  SVEs.  In 
Figure  9,  we  plot  the  cumulative  extreme  value  distributions  of  the  Fatemi-Socie  FIP  as  estimated  over 
grain  sized  volumes.  In  this  plot,  each  datapoint  represents  a  simulated  SVE  of  IN100.  The  scale  of  the 
vertical  axes  was  constructed  such  that  a  true  Gumbel  distribution  would  appear  as  a  straight  line  in  the 
plot.  This  linearization  procedure  and  method  used  to  fit  the  experimental  data  to  a  Gumbel  distribution  is 
outlined  in  Appendix  A.  The  fits  for  these  sampled  distributions  to  the  Gumbel  distribution  appear  in  each 
case  in  Figure  9  as  lines.  As  can  be  seen  in  Table  5,  these  data  are  fit  by  the  Gumbel  distribution  with  high 
confidence  (i.e.,  R2>0.97),  even  for  a  sample  size  of  just  25  simulated  SVEs;  moreover,  the  fitting 
parameters  do  not  change  significantly  with  increased  sample  sizes. 
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Figure  10  shows  the  distributions  of  the  extreme  value  Fatemi-Socie  FIPs  for  the  different  strain 
magnitudes  tested,  0.3%,  0.5%,  and  0.7%.  Again  this  plot  is  linearized  for  the  Gumbel  distribution  as 
outlined  in  Appendix  A.  Additionally,  in  the  extreme  values  of  the  FIPs  are  plotted  in  Figure  10  on  a  log 
scale  in  order  compare  the  distributions  for  the  different  strain  magnitudes  on  the  same  plot.  The 
parameters  for  fit  of  the  Gumbel  distributions  to  these  observations  are  given  in  Table  5.  We  note  that  as 
the  strain  magnitude  increases,  the  extreme  value  distributions  of  the  Fatemi-Socie  FIP  tend  to  move  to  the 
left.  This  correlates  with  the  fact  that  as  the  strain  magnitude  increases  the  driving  forces  for  fatigue 
damage  formation  are  higher.  Thus,  fatigue  cracks  form  more  rapidly  and  the  portion  of  fatigue  life 
attributed  to  fatigue  crack  formation  is  shorter. 

Examination  of  the  parameters  for  the  Gumbel  distribution  fits  of  the  estimated  extreme  value  FIP 
distributions  are  informative  in  terms  of  the  relative  dispersion  of  the  distributions  for  the  different  strain 
magnitudes.  For  the  Gumbel  distribution,  OCn  is  an  inverse  measure  of  the  dispersion  of  the  largest  values 

of  the  initial  population  (i.e.,  the  parent  distribution  of  the  extreme  value  distribution)  while  U  can  be 

considered  as  the  characteristic  largest  value  of  the  initial  population.  Thus,  the  quantity  jlln  is  an 

indication  of  the  dispersion  of  the  extreme  value  distribution  relative  to  the  characteristic  largest  values  of 
the  initial  population.  In  other  words,  larger  values  of  jl ln  indicate  greater  variation  around  the 

characteristic  largest  values  than  smaller  values.  In  our  case,  as  the  strain  magnitude  increases  the  variation 
of  the  extreme  value  distributions  of  the  distributions  (i.e.,  jun  from  Table  5)  tends  to  decrease. 

These  data  corroborate  the  S-N  curve  for  IN100,  as  seen  in  Figure  2(a-b)  with  the  noted  exception  at 
1150MPa  where  there  are  arguably  two  competing  modes  of  failure.  As  the  stress/strain  magnitude 
increases  the  overall  variation  in  the  fatigue  lives  decrease.  Our  results  similarly  show  that  variation  of 
fatigue  life  associated  with  crack  formation  decreases  as  the  magnitude  of  applied  cyclic  straining 
increases. 

To  understand  how  local  attributes  of  microstructure  influence  the  extreme  value  fatigue  response  as 
estimated  via  the  FIPs,  we  consider  the  microstructure  attributes  in  the  neighborhoods  where  the  extreme 
value  FIPs  were  recorded.  Here  we  consider  two  forms  of  the  Schmid  factor.  The  apparent  Schmid  factor 
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for  orientation  g,  is  based  purely  on  geometry,  i.e.,  =  max{cos(^')cos(/l'),!  =\...n\  >  where  ft  is  the 

angle  between  the  ith  slip  plane  normal  and  loading  direction  and  X  is  the  angle  between  the  ith  slip 
direction  and  the  direction  of  uniaxial  applied  stress  on  the  SVE.  The  total  number  of  slip  systems  is  n.  The 
local  Schmid  factor  calculated  over  a  volume  to,  m? ,  is  computed  as  the  ratio  between  the  maximum 
resolved  shear  stress  (as  estimated  from  the  FE  simulations)  divided  by  the  global  SVE  uniaxial  applied 
stress,  i.e.,  mf  =  max ( t‘rss j j a  for  i=l...n.  In  contrast  to  the  apparent  Schmid  factor,  the  local  Schmid 

factor  considers  the  local  stress  state  as  affected  by  intergranular  interactions  and  microplasticity.  For  the 
calculations  given  below,  the  volumes  over  which  we  calculate  the  local  Schmid  factor  correspond  to  the 
same  volumes  over  which  we  calculate  the  FIPs.  In  contrast,  the  apparent  Schmid  factor  is  based  on  the 
orientation  and  relative  loading  direction  at  a  single  location  (i.e.,  for  a  particular  element)  or  for  a  grain. 

First,  we  focus  directly  on  the  zones  where  the  extreme  value  FIPs  has  been  identified.  In  Figure  11, 
we  plot  the  Fatemi-Socie  FIP  as  estimated  over  the  volume  of  a  single  element  versus  the  corresponding 
local  Schmid  factor  for  that  same  volume  for  the  applied  maximum  strain  magnitudes  of  0.3%,  0.5%  and 
0.7%,  respectively.  It  is  noted  that  as  the  maximum  cyclic  strain  amplitude  increases,  the  local  Schmid 
factors  associated  with  the  extreme  value  FIPs  tend  to  decrease  slightly.  This  is  logical  because  even 
though  the  far  field  stress  may  be  increasing  the  flow  stress  will  not  change  significantly.  Thus,  the  local 
Schmid  factor  will  decrease  in  magnitude  as  the  far  field  load  increases.  Secondly,  we  point  out  that  the 
extreme  value  FIPs  seen  in  Figure  1 1  are  primarily  associated  with  cube  slip.  Note  that  in  Figure  1 1  we  plot 
both  the  local  Schmid  factor  for  octahedral  slip  and  cube  slip  in  each  grain  with  the  extreme  value  FIP  for 
each  SVE  simulated  for  0.3%,  0.5%,  and  0.7%  strain.  In  other  words,  the  grains  in  which  the  extreme  value 
FIPs  are  identified  all  appear  to  be  unfavorably  oriented  for  octahedral  slip. 

We  also  consider  certain  correlated  microstructure  attributes  at  the  locations  of  the  FIPs  as  the  second 
part  of  the  microstructure-sensitive  extreme  value  statistical  framework.  As  discussed  previously,  the 
primary  microstructure  attribute  considered  in  this  study  is  crystallographic  orientation.  We  construct  the 
extreme  value  marked  radial  correlation  function  for  the  apparent  Schmid  Factor  (i.e., 
(mga  ,m*  |r,Q))-  By  using  the  apparent  Schmid  factor,  we  really  are  directly  quantifying  the 

probabilities  of  orientations  and  misorientations  of  slip  systems.  As  will  be  seen  later,  constructing  the 
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extreme  value  correlation  function  with  the  apparent  Schmid  factor  instead  of  the  more  direct 
microstructure  attribute  of  crystalline  orientation  simplifies  the  interpretation  of  the  results.  We  choose  not 
to  use  the  local  Schmid  factor  computed  using  the  pointwise  stress  tensor  acting  on  slip  systems  at  each 
point  in  the  mesh  because  the  local  Schmid  factor  is  more  of  a  response  parameter  than  a  direct  description 
of  the  local  microstructure  attribute  of  interest  (i.e.,  grain  orientation). 

In  Figure  12,  we  plot  the  complete  radial  distribution  estimated  over  the  entire  ensemble  of  simulated 
SVEs  and  the  extreme  value  marked  distribution  functions  sampled  only  at  the  location  of  the  extreme 
value  FIPs  for  the  100  SVEs  cycled  at  0.5%  strain.  The  particular  correlations  plotted  here  for  the  local 
Schmid  factors  for  mg  =0.45  to  0.5  (for  cube  slip)  and  ms  = 0-45  to  0.5  (for  cube  slip)  in  Figure  12(a)  and 

complete  radial  and  extreme  value  marked  radial  distribution  of  Schmid  factors  mg=  0.45  to  0.5  (for  cube 
slip)  and  »ij'=0.45  to  0.5  (for  octahedral  slip)  in  Figure  12(b)  were  selected  because  they  occurred  with  the 
highest  probability  at  the  extreme  value  location  of  the  FIPs.  In  Figure  12(a),  we  look  at  the  auto¬ 
correlations  between  Schmid  Factors  mg=  0.45  to  0.5  (for  cube  slip).  We  see  that  below  one  average  grain 

size,  the  probability  of  finding  mg=0A5  to  0.5  (for  cube  slip)  at  the  location  of  the  extreme  value  response 
is  as  high  as  33%;  in  contrast,  the  probability  of  finding  this  particular  orientation  in  the  complete 
microstructure  ensemble  is  only  20%.  Similarly,  the  probability  of  finding  msa=  0.45  to  0.5  (for  cube  slip) 

near  wf  =0.45  to  0.5  (for  octahedral  slip)  (as  can  be  seen  in  Figure  12(b))  is  as  high  as  0.47  at  a  distance  of 

one  average  grain  size;  whereas,  the  probability  of  these  correlated  orientation  existing  in  the  complete 
microstructure  ensemble  is  less  than  0.2.  In  essence,  we  observe  a  high  probability  of  finding  a  high  FIP  in 
a  region  with  clusters  of  grains  oriented  for  cube  slip  or  clusters  of  grains  oriented  for  cube  slip  surrounded 
by  soft  grains  oriented  for  easy  octahedral  slip. 

With  an  understanding  of  how  the  slip  planes  are  oriented  relative  to  the  loading  direction  as  described 
by  the  apparent  Schmid  factor,  it  is  also  interesting  to  explore  the  distributions  of  the  local  Schmid  factors 
in  the  regions  identified  to  be  critical  by  the  extreme  value  FIPs.  In  Figure  13,  we  plot  the  probability 
density  of  the  local  Schmid  factors  nil  for  all  the  grains  and  for  the  grains  with  the  apparent  Schmid  factors 
ma  =  0.45  to  ma  =  0.5  for  both  the  octahedral  and  cube  slip  systems,  respectively.  When  we  consider  the 
distribution  of  the  local  Schmid  factor  over  all  of  the  grains  we  notice  that  the  maximum  local  Schmid 
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factors  for  octahedral  slip  are  less  than  the  maximum  local  Schmid  factors  for  cube  slip.  This  suggests  that 
the  resolved  shear  stresses  on  the  cube  slip  systems  can  be  much  larger  than  they  are  on  the  octahedral  slip 
systems.  As  expected,  when  the  apparent  Schmid  factor  is  larger  (i.e.,  between  0.45  and  0.5)  for  cube  slip, 
the  corresponding  local  Schmid  factors  are  much  more  likely  to  be  larger  than  the  local  Schmid  factors  for 
cube  slip.  It  is  interesting  to  note,  however,  that  for  larger  Schmid  factors  for  octahedral  slip  (i.e.,  0.45  to 
0.5)  we  observe  that  the  local  Schmid  factor  for  cube  slip  can  be  larger  than  the  local  Schmid  factor  for 
octahedral  slip.  Thus,  in  some  cases  even  when  the  Schmid  factors  for  octahedral  slip  are  high,  the  actual 
resolved  stresses  can  be  larger  on  some  of  the  cube  slip  planes  than  on  the  octahedral  slip  planes. 

4.  Discussion 

The  extreme  value  distributions  of  the  Fatemi-Socie  FIP  are  well  fit  by  the  Gumbel  distribution  (i.e., 
R2>0.97)  in  the  Ni-base  superalloy,  1N100,  in  F1CF.  This  observation  is  evident  with  as  few  as  25 
simulation  points  as  seen  in  Figure  9.  We  hypothesize  that  the  distributions  of  the  extreme  value  Fatemi- 
Socie  parameter  is  directly  related  to  distributions  of  crystallographic  fatigue  damage  formation  in  this  and 
similar  material  systems.  Thus,  assuming  an  averaging  volume  for  the  relevant  FIP  that  corresponds  to 

crack  incubation  at  this  scale,  one  may  postulate  a  relation  such  as  PFS  =  y  .  (2Nmc)'  ,  where  f  f  is  a  fatigue 
ductility  coefficient  appropriate  to  crack  formation  at  the  scale  considered.  This  is  obviously  related  to  the 
form  of  the  Coffin-Manson  law  as  defined  in  Equation  1,  with  rescaling  of  parameter  y  t  to  correspond  to 

a  crack  size  on  the  order  of  microstructure.  Using  this  type  of  relation,  optimization  problems  can  be 
envisioned  regarding  polycrystalline  orientation  distributions  to  achieve  target  HCF  lives  or  variability 
thereof. 

This  work  has  considered  only  crystallographic  fatigue  damage  formation  and  did  not  account  for  the 
effect  of  voids  or  non-metallic  inclusions  relative  to  their  influence  on  the  local  driving  forces  for  fatigue 
damage  formation.  Indeed,  in  IN  100  the  role  of  non-metallic  inclusions  is  very  important  in  forming  life- 
limiting  fatigue  cracks.  Each  of  these  attributes  can  have  an  effect  on  overall  variability  of  fatigue  damage 
formation  in  this  material  system  and  must  be  considered  for  any  new  material  development.  The 
understanding  of  extreme  value  slip  behavior  in  this  material  system,  however,  is  very  important  to  the 
understanding  of  fatigue  crack  formation  around  inclusions/voids.  For  example,  one  can  imagine  a  sample 
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of  IN100  with  two  inclusions  of  very  similar  character  (i.e.,  phase,  size,  shape)  where  a  fatigue  crack  forms 
at  one  of  the  inclusions  and  not  the  other  when  the  material  is  subjected  to  some  specified  cyclic  loading 
condition.  Likely,  an  understanding  of  the  local  grain  character  (i.e.,  size,  shape,  orientation,  disorientation) 
near  the  inclusion  where  the  fatigue  crack  formed  relative  to  the  inclusion  void  of  fatigue  damage  will  be 
sufficient  to  explain  the  different  responses  in  the  neighborhood  around  the  two  inclusions. 

Consideration  of  the  extreme  value  marked  radial  distributions  of  the  apparent  Schmid  factor  allows  us 
to  readily  observe  the  character  of  slip  that  is  predicted  to  exist  at  the  location  of  the  extreme  value 
response  (i.e.,  FIPs).  Here  we  observed  that  grains  oriented  unfavorably  for  octahedral  slip  or  favorably  for 
cube  slip  are  predicted  to  exist  with  high  probability  at  the  locations  of  extreme  value  fatigue  response  in 
IN  100.  Moreover,  clusters  of  grain  oriented  for  cube  slip  or  clusters  of  grains  oriented  for  cube  [32]  slip 
surrounded  by  grains  oriented  favorable  for  octahedral  slip  are  predicted  to  exist  with  high  probability  at 
the  location  of  the  extreme  value  fatigue  response.  These  simulations  support  the  observations  made  in 
1N100  by  Li  et  al.  [22]  who  observed  fatigue  damage  formation  along  {100}  planes  in  grains  oriented 
unfavorable  for  octahedral  slip.  Thus,  following  our  previous  logic,  we  expect  that  clusters  of  grains 
oriented  unfavorably  for  octahedral  slip  near  inclusions  or  voids  are  more  susceptible  to  fatigue  crack 
formation  than  inclusions  or  voids  that  are  not  located  in  such  clusters. 

As  improved,  next  generation  microstructure-sensitive  constitutive  models  become  available,  they  can 
serve  the  needs  of  the  framework  developed  here.  In  fact,  an  interesting  question  is  the  extent  of  accuracy 
required  of  these  models  to  deliver  useful  information  for  minimum  life  (low  probability  of  failure)  design 
of  the  microstructure  against  fatigue.  It  is  an  open  issue. 

6.  Conclusions 

The  proposed  statistical  framework  that  captures  the  extreme  value  distributions  of  the  FIP  responses 
coupled  with  the  extreme  value  marked  correlation  functions  has  allowed  us  to  assess  the  distribution  of 
driving  force(s)  for  fatigue  damage  formation  (i.e.,  crack  nucleation  and  microstructurally  small  crack 
growth)  in  polycrystalline  IN  100.  In  these  simulations,  we  account  for  the  material  microstructure,  applied 
strain  amplitude  and  strain  state,  and  critically  stressed  volume  of  interest  in  HCF  scenarios.  Specifically, 
we  found  that: 
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•  The  extreme  value  distributions  of  the  Fatemi-Socie  fatigue  indicator  parameters  (FIPs)  that 
estimate  the  driving  forces  for  fatigue  damage  formation  appear  are  fit  with  high  confidence  (i.e., 
R2  >  0.97)  by  the  Gumbel  distribution. 

•  The  shape  of  the  extreme  value  distributions  of  the  Fatemi-Socie  FIPs  appears  to  be  well  defined 
with  as  few  as  25  SVEs  (i.e.,  microstructure  instantiations)  as  demonstrated  by  the  similarity  in  the 
parameters  of  the  fits  of  the  observations  to  the  Gumbel  distributions  for  sample  sizes  of  25,  50,  75 
and  100  SVEs. 

•  These  simulations  predict  that  cube  slip  may  play  an  important  role  in  fatigue  damage  formation 
particularly  when  there  are  multiple  grains  oriented  for  cube  slip  clustered  in  the  same  region  or 
clusters  of  grains  oriented  for  cube  slip  surrounded  by  other  grains  oriented  favorable  for 
octahedral  slip.  These  observations  support  the  results  of  previous  experiments  of  fatigue  in  this 
material  system  [22].  Thus,  as  fatigue  crack  formation  is  dominated  by  the  presence  of  non- 
metallic  inclusions  or  voids  for  this  particular  material  system,  we  expect  that  given  two  inclusions 
of  similar  character  that  fatigue  cracks  will  form  preferentially  near  the  inclusion  that  is 
surrounded  by  grains  or  grain  clusters  unfavorably  oriented  for  octahedral  slip. 
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Appendix  A:  Linearization  and  Fit  of  the  Gumbel  Distribution 

The  cumulative  distribution  function  (CDF)  of  the  Gumbel  [27,  28]  (extreme  value)  distribution  for  the 

largest  values  F  can  be  expressed  as 

Fy„  (yn  )  =  exP  -ea'[y^Un)  (A.l) 

where  Un  is  the  characteristic  largest  value  of  the  initial  variable  X ,  and  OCn  is  an  inverse  measure  of 
dispersion  of  the  largest  value  of  X  .  Flere  n  refers  to  the  size  of  the  samples  of  the  initial  variable  X  in 
the  set  of  distributions  of  X  from  which  Y  is  sampled.  The  expression  in  equation  A.  1  is  the  asymptotic 

form  of  the  Type  I  extreme  value  distribution  of  the  largest  values  and  therefore  is  valid  as  the  sample  size 
fl  approaches  infinity  [28]. 

To  linearize  the  Gumbel  distribution,  we  solve  equation  A.l  for  the  estimated  parameter  y  given  a 
probability  p  ,  i.e., 


(A.2) 


Thus,  we  now  have  an  expression  in  the  form  of  y  =  mx  +  b  where  V  =  I  n  ^  I  n  ( 1  /  /; )  j  ,  m  =  CCn, 
X=y„,  and  b  =  -anUn  . 

Ranking  the  data  from  a  sample  of  the  extreme  value  distribution  K  in  the  order  of  smallest  to  largest 
(i.e.,  y\ , y] y"  ),  one  can  determine  estimators  F^  ( y\  ) ,  Fy*  [y]  'j  ,  F},  [y\  j ,  etc.  using  the  median 
rank  statistic  defined  as 

_  j  -  0.3 

V  =  ~ -  (A. 3) 

'  h  +  0.4 

•  til  •  •  •  ,  . 

for  the  J  ranked  observation  in  a  sample  of  size  ll  [58].  Flere  the  superscripts  on  yn  refer  to  their  rank 
for  a  particular  sample  of  Yn  .  The  procedure  then  to  plot  a  given  dataset  on  a  linearized  scale  for  the 


Gumbel  distribution  is: 


1.)  Using  Equation  A. 3,  we  estimate  the  probability  p  for  each  observation  yn  in  the  dataset. 
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2. )  Determine  In  In  (1/ p  jj  for  each  observation  and  plot  it  against  the  original  observation. 

3. )  Using  a  least  squares  fit,  a  line  of  the  form  y  =  mx  +  b  is  fit  to  the  data. 

4. )  The  estimated  parameters  for  the  distribution  is  determined  according  to  m  =  C/.n  ,  X  =  V  ,  and 

b  =  -anun . 


This  general  procedure  is  outlined  and  explained  in  detail  by  Kapur  and  Lamerson  [58]  for  similar 
distributions. 
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Tables 

Table  1  The  slip  character  of  several  different  single  crystal  Ni-base  superalloys. 


Material 

Volume 
fraction  y' 

Average 
Size  of  y' 
(nm) 

y'  Shape 

Temp.  (C) 

Strain  Rate 

(s'1) 

Slip  Character 

Ref. 

Rene  N4 

0.65 

250 

Cuboidal 

RT,  650, 
760,  870, 
980 

1  to  2. 5x1  O'5 

[111]  oriented  crystals  exhibited  {100}  slip 
traces  at  all  temperatures  and  range  of 
orientations  exhibiting  {100}  slip  traces 
(instead  of  {111})  increased  with 
temperature 

[19] 

NIMONIC 

105 

0.47 

17.2 

Spherical 

10  to  877 

7x1 0'5, 
1.4x1 0'4 

No  observed  {100}  slip  traces,  only  {111} 
type 

[21] 

SC16 

0.35,  0.05 

450,  80 

Cuboidal 

650,  750, 
850 

10'5 

Observed  slip  traces  on  {100}  planes  at 
650C  and  750C  for  <1 1 1  >  oriented  crystals 

[18] 

Heat 

Treated  PW 
alloy 

0.51 

500 

Cuboidal 

RT 

10'5 

[110]  oriented  crystals  exhibited  {111}  slip 
traces,  [111]  oriented  crystal  exhibited 
{100}  slip  traces  and  [100]  oriented  crystals 
exhibited  neither 

[20] 

Table  2  The  slip  character  of  several  different  polycrystalline  Ni-base  superalloys.  Because  in  many  cases 
the  volume  fraction  of  the  y'  precipitates  was  not  reported,  we  list  the  combined  weight  %  of  Al+Ti  which 
are  the  limiting  components  for  the  formation  of  the  y'  precipitates. 

*  In  this  material  we  report  the  combined  wt.  %  of  Al+Nb+Ta  which  are  in  this  case  the  limiting 
components  for  the  formation  of  the  y'  precipitates. 

**  Although  this  value  was  not  reported  in  the  paper  it  was  given  by  the  authors  who  wrote  this  paper. 
Values  not  reported  in  the  respective  references  are  listed  as  “nr”  for  not  reported. 


Material 

Grain 

Size 

(pm) 

wt.% 

Al+Ti 

Volume 

fraction 

Y‘ 

Ave. 

Size  of 

Y’  (nm) 

Y'  Shape 

Temp. 

(C) 

Slip  Character 

Ref. 

MARM 

004 

nr 

10.55* 

0.5 

-1000 

Cuboidal 

RT,  600 

Facets  on  {1 1 1}  at  RT  and  facets 
on  {100}  at  600C 

[59] 

Nimonic 

API 

40-50 

7.4 

nr 

500 

Cuboidal 

RT 

Most  cracking  on  {111},  but 
cracking  on  {100}  observed  when 
slip  is  limited  to  single  grains 

[60] 

Udmit  700 

nr 

7.5 

nr 

200 

Cuboidal 

RT,  850 

Crack  growth  observed  only  on 
{100}  planes  at  both  RT  and 
850C 

[61] 

INI  00 

30 

10.2 

0.6** 

-1000 

Cuboidal 

RT,  538 

At  RT  cracks  form  on  {1 1 1}  and 
at  538C  cracks  form  on  {100} 

[22] 

Rene 

88DT 

26 

5.8 

nr 

100-200 

Spherical 

593 

Cracks  form  and  propagate  along 
{111}  planes 

[23] 
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Table  3  Model  parameters  used  for  all  simulated  microstructure  in  the  mesh  quality  study. 


Edge  Length  of  Microstructure  Volume  40  pm 

Average  cube  root  grain  size  1 0  pm 

Number  of  Grains  77 

Ratio  of  Minimum  Spacing  of  Grain  Centers  to  Average  Distance  Between  Grain 

Centers  0.3 


Table  4  A  description  of  the  different  meshes  considered  in  the  mesh  study  for  Microstructure  A. 


Number  of  Elements 
Along  Edge 

Total  Elements 

Element  Volume 
(mm3) 

Ratio  of  Element  Volume 
to  Average  Grain 
Volume 

8 

512 

1.25x10  7 

0.1250 

10 

1000 

6.40x1 0'8 

0.0640 

12 

1728 

3.70x10 8 

0.0370 

14 

2744 

2.33x1 0"8 

0.0233 

16 

4096 

1.56x1 0'8 

0.0156 

18 

5832 

1.1 0x1 0'8 

0.0110 

20 

8000 

8.00x10 9 

0.0080 

22 

10648 

6.01x10'9 

0.0060 

24 

13824 

4.63x1 0'9 

0.0046 

26 

17576 

3.64x1 0'9 

0.0036 

28 

21952 

2.92x1 0'9 

0.0029 

Table  5  The  parameters  for  the  least  squares  fit  of  the  Gumbel  distribution  (i.e., 


f  (>’,)  =  exp 


—e 


-an  ( yn  -un ) 


)  for  the  extreme  value  distributions  of  the  Fatemi-Socie  FIPs  sampled 


over  25,  50,  75,  and  100  SVEs,  respectively,  for  0.5%  maximum  strain.  These  extreme  value  distributions 
are  plotted  in  Figure  9.  The  goodness  of  fit  is  also  given  in  terms  of  the  R2  value  that  ranges  between  0  and 
1  with  1  indicating  a  perfect  fit. 


#SVEs 

u„ 

R2 

25 

8186.6 

4.375E-04 

0.973 

50 

8345.4 

4.453E-04 

0.988 

75 

8581.8 

4.402E-04 

0.987 

100 

8727.9 

4.307E-04 

0.991 
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Table  6  The  parameters  for  the  least  squares  fit  of  the  Gumbel  distribution  (i.e., 


Fy„  (>’„ )  =  exp 


-e 


n(. Vn~un) 


)  for  the  extreme  value  distributions  of  the  Fatemi-Socie  FIP  for  the 


different  maximum  strain  levels  simulated.  The  goodness  of  fit  is  also  given  in  terms  of  the  R2  value  that 
ranges  between  0  and  1,  with  1  indicating  a  perfect  fit. 


%  Strain 

un  1 

R2 

0.30% 

6.0  xlO10 

4.619  xlO11 

0.361 

0.952 

0.50% 

8727.9 

4.307  xlO-4 

0.266 

0.991 

0.70% 

3478.5 

3.305  xlO-3 

0.087 

0.981 
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Figures 


Figure  1  Schematic  of  zig-zag  {111}  slip  in  the  {100}  channels.  The  thick  arrows  indicate  the  directions  of 
the  dislocation  segments.  The  [001]  direction  is  the  intersection  line  of  the  {111}  planes  and  is  parallel  to 

the  Burgers  vector  of  the  dislocation  loop  b  =  ^[OOl]  (Figure  from  [18]). 
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Figure  2  The  fatigue  variability  of  IN100  from  physical  fatigue  testing  at  various  stress  magnitudes;  a)  The 
mean  versus  life-limiting  behavior  and  (b)  Cumulative  distribution  functions  for  tests  at  different  stress 
magnitudes  (Figure  from  [24]). 
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Figure  3:  SVEs  generated  via  the  Voronoi  tessellation  based  microstructure  generator  for  1N100  with 
associated  FE  mesh.  Note  that  the  meshes  are  voxellated  and  do  not  exactly  correspond  to  the  boundaries  as 
defined  by  the  Voronoi  tessellation. 
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Figure  4:  The  target  and  optimized  grain  size  distribution  of  a  SVE  generated  for  IN  100  where  V  is  the 
actual  grain  volume  and  <V>  is  the  target  average  volume  or  8.0xl0'6mm3  or  a  cube  root  grain  size  of 
~0.02mm. 
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(a) 


(b) 

Figure  5:  The  mesh  and  contours  of  accumulated  plastic  strain  after  three  cycles  at  0.5%  strain  with  R=-l 
for  the  77  grain  model  for  Microstructure  C  with  (a)  14  and  (b)  28  elements  along  each  edge.  The  elements 
belonging  to  the  grain  with  the  extreme  value  FIP  are  highlighted  in  red. 
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-3 


(a) 


-4 


(b) 

Figure  6  Convergence  of  the  extreme  value  Fatemi-Socie  FIP  for  (a)  the  FIP  calculated  over  a  single 
element  and  (b)  the  FIP  calculated  over  all  the  elements  in  a  single  grain.  Although  convergence  was  not 
achieved  over  the  element  sized  averaging  volume,  variation  was  minimal  over  most  of  the  meshes 
simulated  for  grain  sized  averaging  volumes. 
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-3 


Figure  7  Comparison  of  the  grain  averaged  Fatemi-Socie  FIP  of  the  center  grain  oriented  with  its  {0  0  1}, 
{111}  and  {2  1  25}  planes  perpendicular  to  the  loading  directions,  respectively,  versus  the  number  grains 
along  the  edge  of  the  SVE.  These  SVEs  consisted  of  arrays  of  cuboidal  grains  each  grain  dimensioned  9|t.m 
along  an  edge. 
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Figure  8  Apparent  Schmid  factors  based  on  crystallographic  orientation  for  the  (a)  octahedral  and  (b)  the 
cube  slip  systems. 
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Extreme  Value  FIP 

c) 


Extreme  Value  FIP 

d) 

Figure  9  The  extreme  value  distribution  of  the  Fatemi-Socie  FIP  as  estimated  over  grain  sized  averaging 
volumes  for  (a)  25,  (b)  50,  (c)  75,  and  (d)  100  SVEs  for  the  simulations  cycled  at  0.5%  strain.  The  extreme 
value  Fatemi-Socie  FIP  was  selected  to  be  the  FIP  with  the  highest  magnitude  out  of  all  the  FIP  calculated 
over  each  grain. 
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Figure  10  The  extreme  value  distribution  of  the  Fatemi-Socie  FIP  as  estimated  over  a  single  element  and 


over  the  entire  grain  for  the  simulations  cycled  at  0.3%,  0.5%,  and  0.7%  strain.  The  extreme  value  Fatemi- 


Socie  FIP  was  selected  to  be  the  FIP  with  the  highest  magnitude  out  of  all  the  FIP  calculated  over  each 
element  or  over  each  grain  in  each  SVE. 
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X  10 

(c) 

Figure  11  Scatter  plot  of  the  extreme  value  Fatemi-Socie  FIP  (abscissa)  versus  the  local  Schmid  factor 
(ordinate)  for  simulations  over  (a)  50  SVEs  at  0.3%  strain,  (b)  100  SVEs  at  0.5%  strain,  and  (c)  50  SVEs  at 
0.7%  strain. 
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(a) 


(b) 

Figure  12  For  IN100  the  (a)  complete  radial  and  extreme  value  marked  radial  distribution  for  apparent 
Schmid  factors  777^=0.45  to  0.5  (cube  slip)  and  mft  =0.45  to  0.5  (cube  slip)  and  (b)  complete  radial  and 

extreme  value  marked  radial  distribution  of  apparent  Schmid  factors  77?^  =0.45  to  0.5  (cube  slip)  and 
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m‘ f  =0.45  to  0.5  (octahedral  slip)  are  plotted  for  the  100  SVEs  subjected  to  0.5%  strain.  The  distance  r  that 
separates  the  two  orientations  is  normalized  against  the  average  cube  root  grain  size  d  (i.e.,  20  pm). 
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Figure  13  Probability  density  of  the  local  Schmid  factors  W/  for  the  octahedral  and  cube  slip  systems  for  all 
grains  and  for  the  grains  with  the  apparent  Schmid  factors  ma  =  0.45  to  rna  =  0.5  for  both  the  octahedral  and 
cube  slip  systems,  respectively. 
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